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This is the second part of the notes to the course on quantum theory of large systems 
of non-relativistic matter taught by J. Frohlich at the 1994 Les Houches summer school. It 
is devoted to a sketchy exposition of some of the beautiful and important, recent results of 
J.Feldman and E.Trubowitz, and J. Feldman, H. Knorrer, D. Lehmann, J. Magnen, V. Ri- 
vasseau and E. Trubowitz. Their results are about a mathematical analysis of non-relativistic 
many-body theory, in particular of the Landau- Fermi liquid and BCS superconductivity, us- 
ing Wilson's renormalization group methods and the techniques of the 4— expansion. While 
their work is ultimately aimed at a complete mathematical control (beyond perturbative 
expansions) of systems of weakly coupled electron gases at positive density and small or zero 
temperature, we can only illustrate some of their ideas within the context of perturbative 
solutions of Wilson-type renormalization group flow equations (we calculate leading-order 
terms in a -^-expansion, where iV is an energy scale) and of one-loop effective potential cal- 
culations of the BCS superconducting ground state. We therefore urge the reader to consult 
the following original articles : 

[1] J. Feldman and E. Trubowitz, Helv. Phys. Acta 63, 156 (1990) 
G. Benfatto and G. Gallavotti, J. Stat. Phys. 59, 541 (1991) 

[2] J. Feldman and E. Trubowitz, Helv. Phys. Acta 64, 214 (1991) 

[3] J. Feldman, J. Magnen, V. Rivasseau and E. Trubowitz, Helv. Phys. Acta 65, 679 



(1992) 
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[4] J. Feldman, J. Magnen, V. Rivasseau and E. Trubowitz, Europhysics Letters 24, 437 
(1993) 

[5] J. Feldman, J. Magnen, V. Rivasseau and E. Trubowitz, Europhysics Letters 24, 521 
(1993) 

[6] J. Feldman, J. Magnen, V. Rivasseau and E. Trubowitz, Fermionic Many-Body Models, 
in Mathematical Quantum Theory I: Field Theory and Many-Body Theory, J. Feldman, 
R. Froese and L. Rosen eds, CRM Proceedings and Lecture Notes, 1994. 

[7] J. Feldman, D. Lehmann, H. Knorrer and E. Trubowitz, Fermi Liquids in Two Space 
Dimensions, in "Constructive Physics", V. Rivasseau (ed.), Lecture Notes in Physics 
vol. 446, Berlin, Heidelberg, New York: Springer- Verlag, 1995. 

In Chapter 1, we review some of the standard material on thermodynamics, quantum 
statistical mechanics and functional integration. Sources for this material can be found in 

[8] D. Ruelle, Statistical Mechanics (Rigorous results), London and Amsterdam: Benjamin, 
1969. 

[9] O. Bratteli and D. W. Robinson, Operator Algebras and Quantum Statistical Mechan- 
ics, Berlin, Heidelberg and New York: Springer- Verlag, 1987 (I, second edition), and 
1981 (II). 

[10] J. W. Negele and H. Orland, Quantum Many-Particle Systems, Frontiers in Physics vol 
68, New York: Addison- Wesley, 1987. 

Refs. 6 and 7 above. 

In Chapter 2, we consider weakly coupled electron gases. We recall the notion of scaling 
limit (large distance scales and low frequencies). We show that, in the scaling limit, systems 
of free electrons at positive density and small temperature in d — 1, 2, 3, ... space dimensions 
can be mapped onto systems of multi-flavour, free chiral Dirac fermions in 1 + 1 space-time 
dimensions, the flavour index corresponding to a direction on the Fermi sphere. This is the 
first manifestation of a "principle of dimensional reduction" which says that, in the scaling 
limit (infrared domain), electron gases at positive density and very small temperatures have 
properties analogous to those of certain 1 + 1 dimensional (relativistic) models of Dirac 
fermions with four-fermion interactions. 

We then introduce the basic ideas underlying a renormalization group analysis of weakly 
coupled electron gases and show that there is a parameter, related to a (dimensionless) 
ratio of energy scales, that plays a role analogous to the number, N, of fermion flavours 
in 1 + 1 dimensional models of Dirac fermions with four-fermion interactions, such as the 
Gross-Neveu model. This enables us to study the renormalization flow in weakly coupled 
electron gases by using -^-techniques. We follow ideas first presented in refs. [1], [4] and [5]. 

Useful additional references to Chapter 2 are: 
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[11] K. G. Wilson and J. B. Kogut, Phys. Reports 12, 7 (1974) 

[12] J. Frohlich, R. Gotschmann and P. A. Marchetti, "The Effective Gauge Field Action of 
a System of Non-Relativistic Electrons", to appear in Commun. Math. Phys. 1995 

(See also Chapter 5 of Part I.) 

In Chapter 3, we derive the renormalization group flow equations, to leading order in a 
expansion, for systems of non-relativistic electrons at positive density and zero temperature 
interacting through weak two-body potentials of short range ("screening"), following ideas 
in refs. [1] - [5]. Special attention is given to understanding the striking similarities with 
1 + 1 dimensional Gross-Neveu type models and to an analysis of the flow of BCS couplings 
("BCS instability"). Besides refs. [1] - [7] and [11], the following references will prove useful: 

[13] A. L. Fetter and J.D. Walecka, Quantum Theory of Many-Particle Systems, New York: 
McGraw-Hill, 1971. 

[14] J. Solyom, Adv. Phys. 28, 201 (1979); (RG methods for one-dimensional systems). 
[15] J. Polchinski, TASI Lectures 1992; (effective actions). 

[16] R. Shankar, Rev. Mod. Phys. 66, 129 (1994); (RG methods for (d > 2) -dimensional, 
non-relativistic many-particle systems). 

[17] D. J. Gross and A. Neveu, Phys. Rev. D 10, 3235 (1974); (^-expansion in the Mitter- 
Weisz- Gross-Neveu model). 

[18] K. Gaw§dzki and A. Kupiainen, Commun. Math. Phys. 102, 1 (1985); (rigorous 
renormalization group analysis of the Gross-Neveu model). 

[19] G. Gallavotti, Rev. Mod. Phys. 57, 471 (1985); (RG methods in scalar field theories). 

[20] J. Polchinski, Nucl. Phys. B 231, 269 (1984) 

In Chapter 4, we study the BCS superconducting ground state, the spontaneous breaking 
of the U(l) global gauge symmetry (the particle number symmetry) in the superconducting 
ground state and the emergence of a massless Goldstone boson associated with the broken 
symmetry. We make use of the Nambu-Gorkov formalism and large- N techniques, (refs. [2] 
- [5]). We emphasize the striking analogies between weakly coupled BCS s-wave supercon- 
ductors and the 1 + 1 dimensional, large-iV chiral Gross-Neveu model. Our discussion of 
symmetry breaking is based on the loop expansion of the effective action of a charged, scalar 
field describing Cooper pairs. In the mean-field approximation, it reduces to the calculation 
of an effective potential that proceeds along the lines of the calculations in ref. [17]. In order 
to understand the dependence of symmetry breaking on dimension and temperature and the 
dynamics of Goldstone bosons, we calculate leading corrections to mean-field theory. 

Useful additional references for Chapter 4 are: 
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[21] P. -G. de Gennes, Superconductivity of Metals and Alloys, New York: Benjamin, 1966. 

[22] J. R. Schrieffer, Theory of Superconductivity, Menlo Park: Benjamin-Cummings, 1964. 

[23] A. A. Abrikosov, L. P. Gorkov and I. E. Dzyaloshinsky, Methods of Quantum Field 
Theory in Statistical Physics, New York: Dover, 1975. 

[24] S. Coleman and E. Weinberg, Phys. Rev. D 7, 1888 (1973); and ref. [17]. 

[25] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961) 

[26] J. Goldstone, Nuovo Cimento 19, 154 (1961); 

J. Goldstone, A. Salam and S. Weinberg, Phys. Rev. 127, 965 (1961) 

[27] D. Mermin and H. Wagner, Phys. Rev. Letters 17, 1133 (1966); 
D. Mermin, J. Math. Phys. 8, 1061 (1967) 

For rigorous results on continuous symmetry breaking, see also 

[28] J. Frohlich, B. Simon and T. Spencer, Commun. Math. Phys. 50, 79 (1976); 

J. Frohlich and T. Spencer, in "Scaling and Self-Similarity in Physics" , 
Progress in Physics, Basel, Boston: Birkhauser, 1983. 

Acknowledgements : J. F. thanks J. Magnen and E. Trubowitz for very helpful dis- 
cussions on the results in refs. 1-7. We thank A. Alekseev, P. Boschung, J. Hoppe and 
A. Recknagel for their help in the preparation of these notes. Special thanks are due to R. 
Gotschmann with whom we had countless, illuminating discussions on renormalization group 
methods in many-body theory, and who participated in the work underlying these notes. 

1 Background material 

In this chapter we recall some basic definitions of quantum statistical mechanics, introduce 
ensembles of identical particles, fermions and bosons, and express Euclidean correlation 
functions in terms of functional integrals. 
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1.1 Thermodynamics and quantum statistical mechanics 

We first recall the definition of statistical ensemble. We assume that a quantum mecha- 
nical system is confined to a box A C IE 3 in 3-dimensional Euclidean space. Important 
characteristics of the system are its volume |A| and the particle number N = n|A|, where n 
is the particle density. 

The quantum-mechanical state space of an N-particle system in the box A is a separable 
Hilbert space 1~t^\ its dynamics is given by a selfadjoint Hamilton operator acting on 
T~L^\ whose energy spectrum is discrete and bounded from below. 

In statistical mechanics, one considers two standard ensembles: One either fixes the num- 
ber of particles in the system {canonical ensemble), or one permits the number of particles to 
vary and keeps fixed only its mean value (N) (grand canonical ensemble). Here, N denotes 
the particle number operator, and ((•)) is the expectation in some state of the system. 

Quantities of interest in thermodynamics are thermodynamic potentials. One of the most 
important such potential is the free energy F: 

f3F([3,A,N) = -\ogZ(P,A,N) (1.1) 

where the canonical partition function Z is defined by 

Z(/5,A,iV) = Tr <) (e-K* ) (1.2) 

The free energy refers to the canonical ensemble, since the number of particles is fixed. One 
also introduces a free energy per unit volume (per particle) by setting 

f(P,A,N) = ±-F(0,A,N) ( = 1f(/3,A,A0 , resp. ) (1.3) 

To an ordered sequence, < r n < • • • < r 2 < T\ < f3, of imaginary times, and an n-tuple of 
observables (bounded operators on H^), ai, a 2 . . . a n , one associates a temperature- ordered 
Green (correlation) function 

^Tr (e^ H e^ H a 1 e-^-^ H a 2 . . . a n e^ H ) (1.4) 

(where we have omitted the subscripts A and superscripts N). 
The grand canonical partition function is defined by 

E((3, A,//) = Tr nA (e-^ H -^) = e W A ^ (1.5) 

The coefficient /i is called chemical potential. The thermodynamic potential Q is the so called 
grand potential. We denote by H\ = ©Ar7i^ V ' ) the direct sum of Hilbert spaces corresponding 
to different numbers of identical particles (bosons or fermions). The Hamiltonian H is then 
a direct sum of A^-particle Hamiltonians. In each A^-particle subspace, H coincides with 
ifW (we again omit the subscript A). In nonrelativistic quantum theory, we assume that 
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the particle number operator, N, and the Hamiltonian commute; in relativistic theories this 
is usually not the case because of particle creation-annihilation processes. 

The laws of thermodynamics are described by differential thermodynamic relations. In 
the grand canonical ensemble the most important such relation is 

dn(P,A,n) = (N)dfi + S A dT (1.6) 

The differential, d, is taken with respect to /3 — (where is Boltzmann's constant) and 
fj,; the box A is kept fixed. The coefficient in front of dp, is the mean value (in an equilibrium 
inverse temperature (3 and chemical potential p) of the particle number operator 
already encountered in the definition of the grand canonical ensemble. The coefficient of dT 
is the entropy of the system. 

In the study of macroscopically large systems it is useful to consider the thermody- 
namic limit, A /" E 3 , and to pass from extensive quantities (F, Q,N, . . .) to intensive ones 
(f,p, n, . ..). For example, in the grand canonical ensemble one defines the pressure, p, by 
setting 

pGM= lim T^n^A^) (1.7) 

A/E 3 |A| 

(where A /* IE 3 in the sense of van Hove or Fisher). The particle density and the specific 
entropy are then given by 

dp dp 

n= dJI> S= dT (L8) 
Grand canonical correlation functions are defined for systems in a box A in the same way 
as for the canonical ensemble. They satisfy the KMS (Kubo-Martin-Schwinger) condition 
which we now describe: Given a *algebra A of quantum-mechanical observables, we define 
a state on A as a normalized (p(I) = 1), positive (p(a*a) > 0, Va G A), linear functional 
(p(aa + f3b) = ap(a) + fip{b), Va, b G A, a, (3 G <D). An example of a state is an equilibrium 
state at some inverse temperature (5 defined by 

p p {a) = ^ Tr (e-^-^a) (1.9) 

It is commonly called Gibbs state. Introducing time-dependent observables in the Heisenberg 
picture by setting 

a{t) = jt{H-»N) ae -it{H-nN) ^ aeA (1 10 ) 

we can formulate the KMS condition as follows: For arbitrary elements a and b of A, there 
exists a function F^ b (z) analytic in the strip —f3 < Im z < and continuous on its closure 
such that 

*&{*) = PMt) b) 
F^t - i(3) = Pp (b a(t)) (1.11) 
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Using the cyclic property of the trace 

Tr (ab) = Tr (ba) (1.12) 

one can check that the KMS condition is satisfied for the Gibbs state ( |1.9| ). In the formulation 
fll.lip , the KMS condition can be transferred to the thermodynamic limit. Actually, one can 
view the KMS condition as an equation for equilibrium states of the system. Any state pp 
satisfying ( |1 . 1 1| ) is called an equilibrium state at inverse temperature (3. 

1.2 Systems of identical particles 

In this section we consider a gas of identical quantum-mechanical particles, bosons or 
fermions. The second quantization of a system of identical particles is conventionally de- 
scribed by starting with the Fock space 

n = ©^ =0 ft (A ° (i.i3) 

where is the iV-particle Hilbert space, i.e., 

H (N) = ( W (1))®.JV ^ for bosons 

H (N) = ^ for fermions ( L14 ) 

and H {1) is the Hilbert space of pure states of a single particle. The subscripts s and a 
indicate completely symmetric and antisymmetric tensor products, respectively, according 
to whether the particles are bosons or fermions. Furthermore, (J-C^^ s ° = (7i (1) )® a0 = <D|0), 
where |0) is the vacuum vector in H. 

Let £ denote the position and spin of a particle. Creation- and annihilation-operators, 
and ^(C), are defined as operator- valued distributions on Tt. They satisfy the com- 
mutation relations 



where the commutator [A, B}_ 
the anti-commutator [A, B\ + : 
spin). Moreover, ^(^)|0) = 0. 
The operators 



[* # (0,* # (^)] ± = o 

[*(0,**(»7)] ± = 8(t-ri) 



(1.15) 



= AB — BA is chosen for bosons (integer spin particles), and 
AB + BA is chosen for fermions (particles with half-integer 



*(/) : = /de*(o/(o 

**(/) := JdU(ZW(0 (1-16) 
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/ G H {1) , then generate a *algebra A H m, with (*(/))* = **(/), and 



[*#(/), ¥ # fo)] ± = 
[*(/), = (1.17) 

where (f,g) = J d£ /(C)fl'(O denotes the scalar product on T^W. For fermions, \l/*(/) and 
*(/) are bounded operators: ||^(/)|| = ||**(/)|| = /)■ 

Rather than starting from the Hilbert space 7i and a concrete realization of ( p..!5|) by op- 
erators acting on H, one can start from the abstract *algebra A n a) generated by the identity 
and creation- and annihilation operators ( [L.16| ), with relations (|1.17 ). We obtain represen- 



tations of A n m from states via the so-called GNS (Gelfand-Naimark-Segal) construction, 
as follows: Consider a state p on a *algebra A, i.e., p is a normalized (p(l) = 1), positive 
(p(a*a) > 0, Va G A), linear functional (p(aa + (3b) = ap(a) + /3p(b), Va, b G A, a, (3 G <C) 
on A. Divide A by the left-ideal 

K := {a G -4|p(&a) = 0, V6 G A} (1.18) 

Define a Hilbert space 7i as the norm-closure of V := .4/7C (where {[a], [b]) := p(a*b), i.e., 
||[a]|| 2 := p(a*a)) and a representation fl p of .4. via ([a], n p (c) [5]) := p(a*cb). For the cyclic 
vector Q p := [1] G V, one has that p(c) = (f2 p , n p (c)f2 p ), Vc G A 

The usual Fock space (cf. (|1.13| )) is obtained by choosing p = p , where p is uniquely 
determined by the equations po(a^f(f)) = = po(^*(f)b) , Va, b G A H w. They imply 
Wick's theorem, 

mi f : m^n 

Munnm «•<*»= E £ ,n (/j , ffP(() > : m=n (i-w) 



i=l 

where 



permutations P i=l 



e P = {\ D ' f r |™ (L20) 
1 • " " ' , for fermions 




In this example, the Hilbert space TC obtained from the GNS construction is the Fock space, 
and the cyclic vector fl Po is the vacuum vector |0). 
A useful operation is 

a = II **(fi) - : a ■ (1-21) 

i 

the so-called Wick-ordering. Up to a sign, : a : is the same monomial as a, but with all \l/*'s 
written to the left of all \l/'s. As usual, the sign is equal to 1, for bosons, and to (—1)*, for 
fermions, where t denotes the number of transpositions necessary to move all ^*'s to the left 
of all *'s. 

Remark. In order to become familiar with the GNS-construction, one may consider, as 
a trivial example, the *algebra A := M n ((D) of complex n x n matrices: Choosing p(a) : = 
(x,ax) , where a G M n (C), (., .) denotes the usual scalar product on C n , and x is an 
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arbitrary unit vector in (D n , one easily recovers the representation of M n ((D) on (D n . In this 
example, K, = {a | ax = 0} (consisting of all matrices having their n-th column identically 
equal to zero if x is the n-th basis- vector of (D™), and A/K, = (D n has complex dimension n. 
Another example is obtained by considering a density matrix (e.g. a Gibbs state) on M„((D). 
It leads to a reducible representation of M„((D). 

Returning to A n (i), it is often convenient to choose an orthonormal basis {/i n }^Li in 
TiS l \ and to expand annihilation- and creation operators according to 

oo 

i 

oo 

= E^MO (i-22) 



with 

[*n,*^]± = <W (1-23) 

We assume that every basis element h n , n = 1,2,3, is a C°°- vector for the one-particle 
Hamiltonian H^' . A typical second-quantized Hamiltonian has the form 

H = /^*<«(-^*)<« 

+ s/**te)**«2)V«i,f2)*(&)*te)rffi^2 (1.24) 

where Aa is the Laplacian with, e.g., Dirichlet boundary conditions at <9A, and V is a 
two-body potential. One can rewrite H as 

H = T(tf*,tf) + V(tf*,tf) (1.25) 

where T(\l/*,\1/) = Xi m ,n^m^mn*n * s a positive-definite quadratic form in \1> = (\l/ m ) and 
= (^m) corresponding to the kinetic energy of particles (first term on the right side of 

( |1.24|) ), and V involves terms of higher order in the operators \l/ and and is even in 

and gauge-invariant. 

The time-evolution of a G A n m in the Heisenberg picture is given by 

a = i [H - /iN, a] (1.26) 

(the term proportional to u drops out for gauge-invariant observables) implying that 

d L V 

i ^ n = A nm m m + -— - fiV n (1.27) 

" n 

and 

d R V 

-%K = *mAmn + ~ l*K (1-28) 

where repeated indices are to be summed over, and |^ is obtained from V by (anti-) 

commuting any factor of in V to the very left and then omitting it and is obtained 
by (anti-) commuting any factor of in V to the very right and then omitting it. 
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1.3 Functional integrals: Bosons 

The goal of this section is to rewrite the correlation functions (see Eq.(O)) in terms of func- 



tional integrals. For this purpose we replace the field operators by classical (Grassmann) 
integration variables for bosons (fermions) and integrate the exponential of the action func- 
tional over all configurations of these variables. To obtain well-defined quantities, a cutoff 
has to be introduced (the same for bosons and fermions). Referring to the second quantized 
Hamiltonian introduced previously, we define by replacing A and V by 

A [ K] = f A nm iin,m<K 
nm 1 , otherwise ' 

T/Wf^r* \p) { V a ^ ^ n ^ ces — K (i 3Q\ 

[ , otherwise 

respectively. ^From now on, bosons and fermions are treated separately. For an arbitrary 
operator a on Ti we define 

a( T ) := e ^ w -^) ^H-^) 

where r denotes imaginary time. Thus, we obtain r-dependent fields ^f n (r), ty* (r). We 
replace the bosonic field operators by complex variables 

Mt),K(t)-+Mt)Mt) (1.31) 



satisfying ^( T ) = ^n( r )) with r G [0, f3\. In terms of the complete orthonormal system 
{h n }^ =1 in 7ix\ we understand the ip^s as modes in the decomposition 

^(r,0 = J2Mr)h n (0 (1.32) 

Furthermore, we assume that AM > and VW[i/}*,i/>] > 0, V^*,^- 
We define the Euclidean, i.e., imaginary-time, bosonic action by 



S [Kl {?P*,?P} = [ dr 



o 



T,(r n ^n){r) + H^\r{r)Mr))-^N{r) 



1.33) 



The system is always assumed to be confined to a compact box A. 

As the variable r takes values in the interval [0 , /3], we must specify boundary conditions 
for the variables ?/>(r), ?/>*(r) at r = 0,(3. They are determined by the KMS condition for 
correlation functions. 

Let us consider the correlation function of two bosonic operators \&*(t) and ^(O). By 
analyticity in the complex time variable, the KMS condition implies that 

P/9 (**(t)*(0)) = p^m*(r - (3)) (1.34) 
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The imaginary-time correlation function of operators CLl, . . . Cl n IS defined as 

Mn) ■ ■ ■ a n {r n )) ^ := ^Tr ( e -/»(* w -M*> Ta^n) • • ■ a n (r n )) (1.35) 

The imaginary-time ordering operator "T" rearranges the operators in decreasing time order, 
and the imaginary times Ti,...r n are assumed to be non- negative, pairwise different and 
smaller than (3. Using the KMS condition we obtain the identities 

(**(r)*(0)>^ = P/3 (vT (r)*(0)) = M*(0)**(t - ft) 

= (*(0)**(r - /?) W = (**(r - /3)*(0))^ (1.36) 

This simple consideration implies periodic boundary conditions for bosonic variables: 

r(r + P)= r(r) , + P) = ^{r) (1.37) 

Introducing integration measures 

Vi> = J] #m , ^* = II #m (1-38) 

we have the following result. 

Lemma : The cutoff correlation functions can be written in terras of functional integrals 
as follows: 

= lj P^^e-^^^n^CrOnC,^) (1-39) 

where n^rrij < K. The integral on the r.h.s. is supplied with periodic boundary conditions, 
Eq. ( \1.3^ ). For fi < 0, both sides of the equality are well-defined in the limit k — ► oo ; and 
the equality continues to hold in this limit. 

It is easy to verify this lemma at the level of formal power series, or by appealing to Schwinger- 
Dyson equations. A rigorous proof is somewhat more complicated, but the result is standard. 

1.4 Functional integrals: Fermions 

Our discussion for fermions follows the procedure outlined for bosons, but some room will 
be given to understanding the properties of integration over anticommuting variables. 
Just as in the previous section, we replace the operators ^, ty* by symbols 

y m ,K^Tprn,r n (1-40) 
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where ip n and ip^ are independent Grassmann variables: 

{Mm} = {CO = {Mm} = (1-41) 

To integrate over Grassmann variables, we define the Berezin integral by requiring that 

{## iP*} = = (1.42) 

and 

' "#„V»= = ! (1-43) 



# n l = j dip* 1 = (1.44) 

We recall some properties of the Berezin integral that follow easily from the definition and 
standard combinatorics: Introducing V^ip = dip n . . . dipi and V^ip* = dip* . . . dip* n , we have 
that 

1 = n / diPi ^ n / dip* ip* = J V^iP*V [n] iP A... iP n iP* n ■ • • "01 (1-45) 

1=1 j=l 

It is easy to see that, for M^[ip*, ip\ a monomial in ip 1 , . . . , ip n , ip^, . . . , 

V [n] iP*V [n] iP M [n] [ip*, ip} = (1.46) 

unless AfM[V>*, ^] = A-01 • • ■ V'nV'n • • ■ ^, A G C. From Eqs. ( |L45|) and ( |Q6| ) and from the 
standard definition of determinants we derive that 

v [n] ip*v [n] ip e -(r,K^) = det K [n] ( L47 ) 

where is a regular n x n matrix. In fact, this determinant formula may be used as a 
basis for a definition of integration over Grassmann variables. 
It is convenient to introduce the left derivative by setting 

5 -ip h ...ip ik ...ip in = (-l)( fe+1 ) \^—ip ik )ip il ...ip ik _ 1 ipi k+1 ...ip in 



= (-1)^+%...^^...^ (1.48) 

where we assume that ii, . . . i n are pairwise different. The derivative is defined similarly. 
All anticommutators involving derivatives vanish, except for 

{fe 4 }=w 'i }=% (L49) 

These rules yield an integration by parts formula for Berezin integrals. Let F^[ip*,ip] 
and [ip*, ip] be two monomials in ipi, . . . , ip n , ip^, . . . ,ip* of degrees / and g. Then 

(_l)/+i f v^iP*V^iP 6F ^' r] G[iP, r] (1-50) 

J dlp m 

12 



A similar formula holds when is replaced by g^r- In particular, choosing G to be the 
exponential exp{ — (ip*, K^ijj)}, we arrive at the following identity 

(K^-\ m ( pW^ e -W^) ^[^^1 (1.51) 
J dip* 



We define the Euclidean action for fermions by the same formula as in Eq. ( |1.33| ). We 



find the boundary conditions at r = 0, (5 for Grassmann variables by using the KMS 

condition. The calculation proceeds along the same lines as for commuting variables, but 
we get an extra minus sign in the last equality of ( |1.36| ), due to fermionic imaginary-time 



ordering. Thus, we conclude that the boundary conditions for fermionic variables are anti- 
periodic: 

r(r + /3) = -V*(r) , V(r + (3) = -V(r) (1.52) 
With these boundary conditions, formula ( |1.39|) holds true for fermions. The proof is much 



easier than the one for bosons, since, after imposing a cutoff, it reduces to pure multi-linear 
algebra involving ( |1.50| ) and ( |1.51| ). 

In the following chapters we shall not explicitly refer to the finite- volume (infrared) cutoff 
A and the ultraviolet cutoff (k < oo) anymore, and we shall work with functional integrals in 
a formal way. In most instances it is, however, straightforward to justify our manipulations. 
For simplicity, we shall study systems at zero temperature; but our method can be used for 
positive temperatures as well. 

We shall assume that the reader is familiar with basic notions and results in the quantum 
theory of systems of non-relativistic, non-interacting fermions at positive density, such as the 
Fermi sphere and -surface, the Fermi momentum h,F and the Fermi velocity vf = — ■ We 
shall use units such that h = 1 (and thus shall not distinguish between "wave vector" and 
"momentum" ) . 



2 Weakly coupled electron gases 

On a microscopic scale (~ 10~ 10 m), many systems of condensed matter physics can be de- 
scribed approximately in terms of non-relativistic electrons — which are fermions — with 
two-body interactions, moving in a static background. We are interested in studying such 
systems in thermal equilibrium at some temperature T and chemical potential fi. The 
Heisenberg equations of motion and the equations for equilibrium states (KMS condition) 
of the microscopic system are, i.g., not exactly solvable. Our main interest lies, however, in 
predicting transport quantities like conductivity, which only depend on physical properties 
of the system at mesoscopic length scales (~ 10 _6 m) characterized by a dimensionless pa- 
rameter A 3> 1 (to be thought of as the ratio of meso- to microscopic scale). Such quantities 
are therefore calculable from processes involving momenta of order /cf/A around the Fermi 
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surface, i.e., from properties of the scaling limit (large A, low frequencies) of the system. One 
can try to extract information on the scaling limit of the system without explicitly solving 
the microscopic equations. Techniques that sometimes allow one to do this involve a prin- 
ciple of dimensional reduction — i.e., the observation that, in the scaling limit, systems of 
non-relativistic (free) electrons in d spatial dimensions behave like systems of multi-flavour 
Dirac fermions in 1+1 dimensions — as well as renormalization-group-improved perturbation 
theory around the non-interacting electron gas. These techniques will be explained in the 
following chapters. 



2.1 Free electrons and dimensional reduction 

In this section, we show how a system of non-relativistic free electrons in <i-dimensional space 
can be approximated by multi-flavour relativistic fermions in 1+1 dimensional space-time. 
We start from the action 

S (r^)=Y.Jd d+1 xrAx)(d ~^-riMx) (2-1) 

k 2 

with a prescribed chemical potential /x (related to the Fermi momentum kp by /i = and 
to the particle density by n = iT^kp / {2n) d , where is the volume of the d-dimensional unit 
ball, and the factor 2 accounts for the spin orientations). For simplicity, we work at zero 
temperature. 

The Euclidean free fermion Green's function is given by (x = (t, x), y = (s,y), t > s, 
with t and s now denoting imaginary time) 



G° aa ,(x-y) = (e tH V a (x)e-^ H %,(y)e- sH ) 



Sa*> I d d+1 k- -p (2.2) 



e 



-iko (t—s)+ik(x—y) 



ik o - (2^ - 
k 



where 



= s m , ld d ke^-iiy-''-"^-''h' 1 ^ (2.3) 

The equality (|2.3| ) follows with the help of the residue theorem. In the following, we are 
interested in physics at a mesoscopic length scale, or, more precisely, in the scaling limit 
(very low momenta and frequencies). Then it suffices to determine the leading contribution 
to G G aa ,{x — y) at arguments x and y with vp\t — s\ + \x — y\ ~ X/kp (where vp = ^ is 
the Fermi velocity). This contribution comes from modes with momenta in a shell S { p ] of 
thickness kp/X around the Fermi sphere Sp: 



S F = {k e H d \ k 2 = k 2 F } 

A < 

2X 



S { F X) = \k g R d | dist(jfe, S F ) < — 
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States of low energy describe electron-hole pairs with momenta in the vicinity of the Fermi 
surface of the system. 

Let us rewrite ( |2.2p in terms of new variables u, p\\ and p , with k F Q G Sp, Po = k , and 
k = (kp + p\\)u. If k G then \p\\\ <C kp, and we can approximate the integrand of ( j2.2| ) 
by dropping the p| term in the denominator. For large A and small frequencies, 

GJU*-y)«*«*< / ^r^f 1 e ik ^ s ^G c {t - s,u{x- y)) (2.5) 

where da(u) is the uniform measure on the unit sphere, and G c is the two-dimensional, chiral 
propagator 

G c (t, = - / (Jpotfpu (2.6) 

J ip - v F pn 

We will set Vp = 1 in the remainder of this section. Introducing the complex variable 
z = it + £ and its complex conjugate z, it is easy to verify that G c (z) satisfies 

-2idG c {z) = 5 {2) {z) 

In other words, G c is the Green's function of a chiral Dirac fermion in 1+1 dimensions. 
Explicitly, 

G c (z) = — (2.7) 

71 Z 



One can further approximate the cJ-integration in Eq. (|2.5f) by a summation over discrete 
directions ujj by dividing 5^ into iV small boxes B^., j = 1, . . . N, of roughly cubical shape: 
is centered at aJj G Sf and has approximate side length kp/X; note that iV w const A d_1 . 
Thus 

G° aal (x-y) w -5 ff(7 ^e^^ / ip d P]l dp ± : (2.8) 

^ J *Po - P|| 

where p = pyo; + pj_ is a vector in B$. — fc^tJj, and p G H. 

We have shown that, in the scaling limit, the behaviour of a (i- dimensional, non-relativistic, 
non-interacting electron gas is described by iV flavours of free, chiral Dirac fermions in a 1 + 1 
dimensional space-time. The statement is to be understood in an appropriate sense, since 
the propagator G c (t — s, u(x — y)) actually depends on the "flavour index" lo. But the energy 
of an electron or a hole with momentum k only depends on pn, where p\\ = kuo — kp, and 



it is proportional to \p\\ | , just as for relativistic fermions in 1 + 1 dimensions. 



2.2 Weakly coupled electrons and the renormalization group (RG) 

We have presented a somewhat unusual description of the free fermion system, because 
expressions such as (|2.8|) provide a convenient starting point for treating interacting fermions 
by a perturbation expansion. In particular, one may hope that interesting physical quantities 
have an expansion in powers of 1/A analogous to the 1/N expansion in the Gross-Neveu 
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model. Later we will see that this is the case for weakly coupled systems, and, as in the 
Gross-Neveu model, the large number of flavours N encountered in the mesoscopic regime 
will play an important role in the actual calculations. We will, however, also see that 
perturbation expansions usually cannot be applied "naively", but have to be improved by 
applying renormalization group (RG) methods. 

We consider systems with a Euclidean action of the form 

S(P,if>)=S (P,if>) + S I {P,if>) 

where S (ili*, ip) is the quadratic term fl2.1|) , and Si(ip*, ip) = g P(ip*, ip) is some higher order 
polynomial interaction to be specified below; we assume that the dimensionless coupling 
constant go is small. Since correlation functions of the fermion fields do not have 

good scaling behaviour, we first split off the oscillatory factors associated with the "direction 
index" Q (see, e.g., Eq. ( |2.8|) ) by expanding the fermions 

V(x)=Y / e l ^dx) (2-9) 

into quasi-particle operators ty$(x). We are interested in calculating connected correlators 
(temperature-ordered Green functions) of the form 

(A a ) 2 "(^ 1 (Aa; 1 ) • ■ ■ ^(Aa^^'i) • • • ^(A<)>* (2.10) 

at large distance scales A; the factor in front of the bracket, involving the scaling dimension a 
of the quasi-particle fields, accounts for the scaling behaviour of the quasi-particle operators. 
Strictly speaking, instead of considering the expectation values above, one ought to "smear 
out" the fields with some test function (an approximate 5- function) h(£): 

*l» = (A)- ((m) / d d+1 yh((x -fl/A)*f(y) 

Because we are interested in the large-scale behaviour, we can work with test functions h 
of compact support in momentum space, e.g. with supp h = {k e IR^ 1 ) k 2 < k 2 F }. After 
Fourier transformation, the connected correlation functions of smeared quasi-particle fields 
are given by 

(A"-*- 1 )*^) • • • h(k' n ) Glih/X, . . . k'JX) 

with 

Gf n = (^(h) ■ ■ ■ *a n (k n )%,(k[) ■ ■ ■ % n { k 'n))l 
According to Chapter 1, we can express this correlator in terms of a functional integral 

Gl n = | J VrV^e-^^ih) ■ -4hSK) 
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Figure 1: Integrating out degrees of freedom 

where kj G Bq. — k F ujj, etc. (because of the restriction on the support of the test functions). 
We now decompose the fermion fields into "slow" and "fast" modes by writing t/j — *0< + V>>, 
with 

supp^> c R x (IR d \S> A) ) (region >) 
supp ■0< cRx (region <) 

The large-scale behaviour of the system is described by the effective action S e fi($< : ■0<) 
given by 

e -s eff (V^<) = _L J p^*^ >e - 5 ^>^<^>^<) (2.11) 

where we perform the functional integral over fast modes only; the normalization factor S> 
is chosen so as to ensure that <S , e ff(0, 0) = 0. Obviously, S e ff depends on our choice of the 
scale A. 

The quadratic part, So, of the action splits into two pieces: 

sS m ,$) = s 0t >(fe,$>) + s 0t< (p < ,j, < ) 

The first term and the normalization factor can be absorbed into the definition of the inte- 
gration measure 

rfP(^*,^>) = ^^;^ >e - 5o ->^>^>) 

and we have to calculate 

= exp{-S ,< - (Sj) G0> + ^(Sj; S I ) G o > - S r , Sj) Gi - . . .} (2.12) 
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The abbreviations 



(A;B) := (AB) — (A)(B) 
(A;B;C) := (ABC)-(A;B)(C)-(C;A)(B)-(B;C)(A)-(A)(B)(C) 

etc., denote connected correlators; the subscript " indicates that the expectations (( . )) G o 
are calculated with the help of the infrared-cutoff free fermion propagator G> , in accordance 
with the functional measure gLP^!,, ip>). The second equality in ( |2.12| ) states that the ef- 
fective action is given by sums over connected diagrams. This is the so-called linked cluster 



theorem. Eq. (|2.12|) also shows that S c g contains i.g. far more interaction terms than <Sj; 
for weakly coupled systems, however, the original couplings will remain the dominant ones. 
The proofs of equation (|2.12j ) and of the linked cluster theorem are standard. We therefore 



omit them. Feynman rules for the computation of S e s are easily derived (we do not present 
the details). 



In order to analyze the perturbation series in fl2.12|) further, we must specify the inter- 



actions in the system: We assume that the leading term in Si is a two-body interaction 

Si(r,^) = ^k^Y, J d d+1 xd d+1 y : rMM*)<x - - y°WAy)My) ■■ (2.13) 

where v(x — y) is a smooth short-range potential (we choose units such that h — 1; the 
factor kp~ d then ensures that g is dimensionless) . It is useful to estimate how close to Sp 
one can get with "naive" perturbation theory: The free propagator, — (ik — uk)' 1 , is regular 
in region >, it is in fact of order \/kp. The size of the two-body interaction associated with 
the graph 



(2.14) 



is therefore roughly equal to \go\ 2 \ (which is dimensionless). The wavy line represents the 
interaction potential which is of order go, the oriented solid lines represent electron "half- 
propagators" each of which is of order A^. As long as A 1/y^o, there are no convergence 
problems with the perturbation series — except that the number of diagrams seemingly 
grows too fast; but their relative signs ensure appropriate cancellations (Pauli principle). 

This crude estimate shows that we cannot pass to the limit A —>■ oo without meeting 
infrared divergence problems in the perturbation series. The method to control the scaling 
limit is RG improved perturbation theory. It allows us to analyze the large-scale behaviour 
of the system by carefully keeping track of the (relative) growth of the various terms in 
the effective action (quadratic part as well as couplings) when the integration over fermion 
modes approaches those near the Fermi surface. In fact, for the interaction Q2.13| ), the actual 
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Figure 2: The renormalization group procedure 



RG computations are not terribly involved, since S e s has an expansion in powers of 1/A, and 
focusing on the leading order in ~ drastically reduces the number of diagrams that have to 
be calculated. 

We shall apply an iterative RG scheme patterned on Wilson's approach. In the "zeroth" 
step, we choose some large scale Ao such that Ao <C Vv^o an d calculate the corresponding 
effective action perturbatively, as in fl2.12|) — to leading order in 1/Ao- The action 
depends on a collection of modes corresponding to wave vectors that are located in a shell 
of width around the Fermi surface Sf- Although not essential for our method, we divide 
this shell into iV « const Aq _1 sectors (boxes), as in the previous section. Furthermore, we 
rescale all momenta so that, instead of belonging to B^., they are contained in boxes B$. of 
side length « kp- Our RG procedure consists in iterating the following two steps: 

Dec Choose some (fixed) integer M > 1 and integrate over the fermion modes corresponding 
to wave vectors in Sp \ Sp (where A = M J Ao, for some j). 

Resc Divide each sector in 5 , ^ MA - ) into M d ~ x new sectors to restore the cubical shape of the 
sectors. Then rescale all momenta by k i — > k = Mk. 

To determine the RG flow of the couplings, we have to find out how the various terms in 
the action transform under rescaling and which diagrams contribute, during the integration 
process, to leading order in 1/A. 



19 



3 The renormalization group flow 



In this chapter we attempt to explicitly implement the renormalization group procedure 
described in the last chapter. Our aim is to understand the renormalization group flow to 
lowest order in j- and for small values of g . 

3.1 Scaling of action and fields 

We shall begin with the discussion of the second step and of the scaling dimensions of the 
different terms in the effective action. Since the decimation of degrees of freedom is most ex- 
plicitly described in momentum-space, our renormalization group procedure is implemented 
in the momentum-space representation. Thus, we first calculate the action in terms of the 
sector fields discussed in the previous chapter. 

Let A ^> 1, and assume that the degrees of freedom not lying in the shell have already 
been integrated out. We divide into boxes of roughly cubical shape and approximate 
side length ^p. The number of such boxes is N ~ Ua-iX^ 1 , where tOd-i is the surface volume 
of the unit sphere in d spatial dimensions. Let B$, \uj\ = 1, denote the box which is centered 
at kpu. The support of those modes 

$ a (k) = I d d+1 xe i{ - kot - %£) ^j a (x) (3.1) 



that have not been integrated out, yet, is contained in 1R x Sp^. Sector fields are defined as 

^ a ( x ) : =/ d d+1 ke- i{kot -^- kFQ)s) i)Jk) (3.2) 

J~RxB a 

The Fourier transform of ipj^^x) has support in 1R x (Bq — kpCu), and 

Let us first look at the quadratic part of the effective action of the modes localized in the 
shell S^. We temporarily assume that this part has the form ( j2.1| ) (there will be a finite 
wave function renormalization which we are neglecting at this point). Inserting the Fourier 
transformed fields yields 

= -£/ , d d+1 p^ a (p)U Po -v F ^p + o(^))^Ap) (3-3) 

^ JRx(B s -k F u) \ \A J J 

with vp = In the last line, we have put the term quadratic in p into an error term O(-ry). 
The leading part of the inverse propagator in So is of order O(t)- Since we always restrict 
our analysis to the leading order in i, we can omit this error term in what follows. 
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Next, we consider the quartic term ( 2.13 ). We insert the expansion of the electron fields 
in terms of sector fields and get 

Si(ip*,ip) =g k)T d ^ J d d+1 xd d+1 ye- ikF{i3l - a ^e- ikF ^ 2 - a3) y 



t>l,...u>4 
a. a 1 



^kA x ^ky(y) v ^ - v) 5 ( x o - yo)VW'(2/)^4,<r( x ) ( 3 - 4 ) 

In terms of the Fourier transform of the potential 

v (x) = f d d ke ik *v{k) 



we can write 

v(x)5(x°) = I d d+1 ke- i{kox °- k£) v(k) 



Since we assume that v(x) is a short-range potential, v(k) is a smooth function of k. Applying 
Fourier transformation to the sector fields in (]3.4f), we obtain 



a. a 1 



2vr5(j9i i0 - p Afi - k Q )2ix5(p 2fl - P3,o + k ) 

(2n) d 5(p 1 -p 4 -k + k F (uj! - uj 2 ))(2ii) d 5(p2 - p 3 - k + k F (u 2 - Cu 3 )) 

C 1>CT (Pl)C 2l( 7'(P2)v(^)^3, CT '(P3)^4,a(P4) 

where the ^-integrations range over the boxes — k F uJi. Performing the /c-integration, we 
obtain the expression 



9ok l f d \ £ Jd d+1 Pi---d d+1 P427r6(p lfi +p 2 . 



- P3,0 - Pi,0) 



J 1 ,...u> 4 



(2n) d 5(p 1 +P2-P3-P4 + kpiuj! + UJ2-U3- UJ4)) 

V(P1 ~P4 + kp(u)i - ^))C liCT (pi)C 2iCT ,(p 2 )^3^'(P3)^ 4 , CT (p4) 

Considering the argument of the second (^-distribution appearing in the integrand, the part 
containing the lu^s is of order k Fl while the part involving the p^s is only of order k y. 
Therefore, for the argument of this 5-distribution to vanish, the term containing the cJj's must 
be zero (to lowest order in ^). Furthermore, as the Fourier transform v of the interaction is 
a smooth function, its value under the integral sign is well approximated by v(k F (coi — d} 4 )), 
dropping the box momenta. Consequently we get 

Si = gok F ' d \ £ u(M<3i - ^)) / J d+1 Pi ■ ■ ■ 3 d+1 PA 

(27r) d+1 (5(pi +p 2 -p 3 -P4)^ 1i( 7(Pi)^ 2 , CT '(P2)^3,<t'(P3)^4,<t(P4) (3.5) 

+ terms of higher order in — 

A 
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Instead of studying the quartic term of the original (microscopic) action of the electron gas, 
we should actually study the quartic (and higher- degree) terms of the effective action at scale 

with A = Ao 1 and \go^l\ <C 1- Using cluster expansions to integrate out the degrees 
of freedom labeled by momenta k ^ Sp \ one can show that, for weakly coupled systems, 
the quartic term of the effective action at scale ^ still has the form given in eq. ( |3.5| ), except 
that v(kp{uji — uja)) is replaced by a coupling function g(u)i, ..,cJ 4 ) « v{kp{Qi — £4)), with 
uj\ + UJ2 = 003 + (Moreover, terms of degree larger than four in if)*, ip are very small.) 

Next, we determine the scaling dimensions of action and fields. We rescale the fields in 
such a way, that the supports of the Fourier transformed, rescaled sector fields are boxes of 
roughly cubical shape and with side length approximately equal to k F , and such that the 
quadratic part of the action remains unchanged to leading order in j. The first condition 
implies, that we should scale momenta according to p — p/X. In configuration space this 
corresponds to the scaling x = A£. Thus the scaled sector fields are 

Vw(0 = WoAX) 

where the scaling dimension a of the sector fields still has to be determined. The Fourier 
transformed sector fields are 

Lap) = ^-^Aj) 

The support of a lies in = X(B^ — k F uj). Indeed, B$ is a roughly cubical box with ap- 
proximate side length kp. Inserting the scaled, Fourier-transformed fields into the quadratic 
part (|3.3|) of the action yields 

So = -a 2 ^ 1 -^-^ 1 ) y: j d d+i p^ a (p) Q(<po - v f qp) + o (i)) Lap) 

For the free part of the action to have the same form as for the unsealed fields , we have to 
choose a = |. Thus the scaled fields are 

$oAZ) = ^HaAK) (3-6) 
LAP) = A- ( ^ +1) 4v(p/A) (3-7) 
Let us determine the scaling behaviour of a term in the action of the form 

s n = ^ £ [d d+1 Pl ---d d+1 p 2n 

n. _, f~r^ J 

"IT H"n-"n+lT r"2n 

Cl,...(T2n 

w(p U ■ ■ •P2n)V4, (71 (Pl) • ■ •fe n , CT2n (p2n)(27T) d+1 5(pi + p 2n ) (3.8) 
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The function w is assumed to be homogeneous of degree kgR, i.e. 



w(spi, . . . sp 2n ) = s K w{px, . . .p 2n ) 
Expressing the sector fields in terms of the scaled fields thus yields 

s n = l A - ((n - 1)c! - 1+K) y / d d+ % ■ ■ ■ d d+l p 2r 



J!i h"n = "n + l" 

01,...<X2n 



W(px, . . .p 2n )^c3 u a 1 (Pl) ■ ■ ■' l Pc3 2 n,a 2n (P2n)(s^) + $(pl + p2n) (3.9) 

The exponent of A is called the scaling dimension of S n , i.e., the scaling dimension of S n is 
— ((n — l)d — 1 + k). With n = 2 and k = 0, the quartic term ( |3.5| ) turns out to have scaling 
dimension 1 — d. Furthermore, we see that quartic terms of higher degree in the momenta 
have a smaller scaling dimension. The same is true for contributions to the effective action 
that are of higher degree in the fields. 

The effective action in terms of the rescaled fields thus reads 



S eS = Z- 1 J2j ^ d+1 P^ l(J (p)(^0 - VFUJp^^ip) 



+l^i-i z ~ 2 E ^i,...,£ 4 ) I d d+l p l ---d d+1 p, (3.10) 



2 \ d ~ 

<J,(J f 



^ 1 , CT (Pl)^ 2 , ( 7'(P2)^3, (7 '(P3)^ 4 , (T (P4)(27r) + 5{pi +p 2 -p 3 -Pa 

+ terms of higher order in — 

A 



In expression ( |3.10| ) for S e s we have introduced a wave function renormalization constant Z, 



in order to indicate that the quadratic part of the effective action may flow under renormal- 
ization (decimation of degrees of freedom and rescaling). In the next section we shall derive 
renormalization group flow equations for Z, vp and the coupling constants g(uJi, . . . , UJ4). 
These equations will determine the dependence of Z, vp and g(uJi, . . . , U4) on the scale 
parameter A. 

The expression for S e g on the right hand side of ( |3.10| ) shows that the inverse propagator 

of the sector fields ^aaip) 1S diagonal in the sector index to and that it only depends on po 
and p\\ = up (but not on p± = p — (up)uj). These features are an aspect of the principle of 
dimensional reduction from d + 1 to 1 + 1 dimensions. Indeed, we observe a rather striking 
formal similarity between expression ( ^.10| ) and the action of the Gross-Neveu model of 



interacting, relativistic Dirac fields in 1 + 1 space-time dimensions: The sector index uj plays 
the role of the flavour index of the Dirac fields in the Gross-Neveu model; the number of 
distinct sector indices, ~ const A d_1 , corresponds to the number, N, of flavours of fermions 
in the Gross-Neveu model. The coupling constants, X 1 ~ d g(ui, . . . , dJ 4 ), correspond to the 
coupling constant, $W-, of the quartic term in the Gross-Neveu model. (The correspondence 
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between ^f- and X 1 ~ d g{u, —u, lj', —Co') = X 1 ~ d g B cs{^>-^>') is particularly precise, as discussed 
in Chapter 4.) 

A powerful method to analyze the Gross-Neveu model is the jj- expansion. This suggests 
to analyze non-relativistic, interacting electron gases with the help of a j-- expansion (with 
y ~ jj, for d = 2), and this is precisely what we shall do in the remaining sections, following 
beautiful ideas of Feldman, Magnen, Rivasseau and Trubowitz. In the Gross-Neveu model, 
Z and the velocity of light (corresponding to vp) do not flow under renormalization, to 
leading order in This suggests that, for the electron gas, Z and vp do not flow under 
renormalization to leading order in ^; a prediction that will turn out to be correct! 

In the following sections, we shall always work in momentum space and with rescaled 
sector fields. We shall thus omit the "hat" and the "tilde" from the rescaled fields on 
momentum space. Our analysis will be based on the assumption that Ao ^> 1 and that 
all coupling constants g^(ui, . . . ,004) = g(u\, . . . ,u; 4 ) |a=a ^ 1- We shall determine the 
renormalization flow to leading order in j- (sometimes omitting terms that are of leading 
order in but of high order in g( Q \u\, . . . , Q4)). 

3.2 Integrating out modes 

In the last section, we have understood how the different parts of the effective action of 
an interacting electron gas behave under rescaling. Here we turn to the second step of the 
renormalization group method — the decimation of degrees of freedom. Initially, we assume 
the degrees of freedom in IR d \ to be integrated out, as discussed in Sect. 2.2. In the j'th 
step, the degrees of freedom localized in Sp 43 x °)\s(m j have to be eliminated, where M is 
a positive integer > 2. Thus the scaling factor at scale j is Xj = X M j . At scale we are given 
Z(°\ VpP and functions g( Q \ui, . . . , U4) of the unit vectors oui, . . . , cJ 4 , with Q\ + uj 2 = LV3 + U4. 
As announced, we assume that Ao ^ 1 and max^ 1 ^ 2 ^ 3 {g^°\coi, . . . ,0)4)} <C 1. 

/From this point on, we omit the spin indices when no confusion arises. We propose 
to first discuss the possible intersector scattering geometries, as we would like to better 
understand the structure of the coupling constants g^(ui, . . . , U4). How many independent 
g^(u 1 , . . . , tD 4 )'s exist ? For d = 3 we suppose that uj 3 ^ — cJ 4 . On the unit sphere, there 
are const Aq _1 different u;'s. All choices ui, uo 2 with ui + u 2 = co 3 + lu 4 lie on a cone 

containing cu 3 and u; 4 with symmetry axis UJ3+U4. Therefore there are 0(Xq~ 2 ) choices. Only 
when uj 3 = —0J4, there are as const Aq _1 choices. Similarly, in d — 2, there are exactly 
two choices for uji, uj 2 , if ^3 7^ —^4, and A^ ) m const A choices if cu 3 = —0J4. 

Couplings involving incoming states with uj 3 7^ — 0J4 shall be denoted by g^°\co 1 , . . . , u) 4 ) 
Couplings that involve sector indices with u 3 = —C04 or, equivalently, Q\ = —uj 2 shall be 
denoted by <7bcs(^i>^4) (f° r BCS- scattering). The latter will prove to be crucial for the 
understanding of the superconducting state. We observe that, because of rotation invariance, 
the coupling gees is on ly a function of the angle between £u 1 and oj 4 . 

Technically, we will carry out the decimation of degrees of freedom using perturbation the- 
ory. If, under the renormalization group transformation, g^(ui, . . . , cJ 4 ) = g(ui, . . . , cJ 4 )| A=A3 
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Figure 3: Sector momentum conservation in d = 3 



increases with growing j (and decreasing distance to the Fermi surface), the validity of per- 
turbative results eventually breaks down. 

We start explicit calculations with the renormalization of the electron propagator when 
passing from scale Ao to Ai. We know that every interaction squiggle provides a factor 
\k~ d g(°\ui, . . . C04). What are the dominant radiative corrections of the electron propagator? 
The two possible one loop diagrams are 



(3.11) 



U),p 



u>, p 



u), p 



ui,p 



We call them tadpole and turtle graph. We shall estimate their amplitudes. There is one 
interaction squiggle of order and m const Aq _1 choices for the inner particle sectors, 

denoted by uj . Therefore, both these graphs correspond to contributions of order zero in 
It will prove to be useful to introduce a simplified graphical notation. Namely, we replace 
the squiggle by a point, i.e., each vertex is represented by a cross with two incoming and 
two outgoing lines. The new vertex stands for the sum of the two old vertices that yield the 
same new vertex. For example, the diagrams in (pj.ll[) are now represented by one graph, 
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namely 

uj',k 




(3.12) 



LJ,P LU,p 



What about combinations of these first order corrections? 




(3.13) 



For every vertex corresponding to a factor -^t, there are const A loop sectors to 

choose from. Therefore all these graphs correspond to amplitudes of order one. But, as a 
matter of fact, one observes with amazement that all their amplitudes vanish. The reason is 
that, in all these graphs, there is an oriented loop involving two particle lines with the same 
sector label u and the same box momenta and energies. For each such loop, the integration 
over the loop momentum yields a factor 




in the total amplitude. The two poles of the integrand coincide, and the integrand decays 
quadratically at infinity. The residue theorem then tells us that the integral vanishes. 

Furthermore, we know that graphs which are not one particle irreducible (1PI) yield a 
vanishing amplitude. Indeed, a graph of this type has the following form 




(3.14) 



and one has to sum over uj' and to integrate over k with < \u'k\ < k F . The empty circles 
stand for arbitrary subdiagrams with the indicated external legs. Because of momentum 
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conservation, there is a factor Sq^i S(p — k). Hence only 63' = 63 contributes. But, in this 
case, the argument of the ^-distribution never vanishes since \63p\ < ^j. Consequently the 
amplitude is zero. 

The graphs not taken into account so far are of order 0(1/ Xq) or higher. As an example, 
let us consider the graph 




(3.15) 



Clearly, 63% is determined by 63, 63\ and 63 2 - Thus there remain two summations, over 63\ and 
6j 2 . If 63\ is different from — 63, then the number of uVs satisfying the relations |u;+u; 2 — 63±\ = 1 
is of order 0(Aq~ 2 ). Hence the number of summands is at most of order 0(Ao d_3 ). Since 
each of the squiggles gives a factor t^=t, the amplitude corresponding to this graph is of 

A o 

order O(l/A ), as claimed. 

Conclusion: The corrections of lowest order in j- to the electron propagator arise from 
the tadpole and the turtle graph. In order to identify the renormalized quantities more easily, 
we work with a renormalized action in its most general form (dropping terms of higher order 
in j, A = Ao, Ai, A2, ...) 

S eS = -Z- 1 J2 J d d+1 k^Ak)(iko-v F 63k)^Ak) 

UJ,(T 

-xz-'J: J d d+1 kr^(k)6^Ak)+ o(j) 

LU,(T 

+ \~>h Z ~ 2 ^ g(63 1 ,...63,)Jd d+1 k 1 ---d d+1 k 4 (3.16) 



/ 

<J,<7 



^ 1 ,.(^i)C'(^)^ 3 y( fc 3)fe,.(fc4)(27r) d+1 5(A; 1 + k 2 - k 3 - k 4 



+ 0(y d ) (3.17) 

On the right side of ( |3.16| ), the wave vectors k = (ko, k) are constrained to belong to 1R x B$, 
for a given sector momentum 63. The factor A in front of the second term expresses the fact 
that this term has scaling dimension 1. Denoting the inverse propagator by T, the sum of 
all quadratic terms in the effective action can be written as 

X; j d d+i kr^( k ) T dk)^Ak) (3.18) 

where 

r a (0) = -XZ^d/j (3.19) 

■^-T a (p) |p= = -iZ- 1 (3.20) 

63(V P TM \p=o) = v F Z- 1 (3.21) 
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Because the tadpole and the turtle graph correspond to p-independent amplitudes, there are 
no corrections to Z and Vp in leading order in 



Z (o) 
,(°) 



O((^ (0) ) 2 /Ao) 



(3.22) 

JF +O((gWy/\ ) (3.23) 
Let us compute the amplitude of the tadpole graph renormalizing the chemical potential: 



An" 



T 9 (0) (u,u,uj',u) J dk ± J^ 



dk\\ 



■<\ku\<ki 



lim 

t|0 



dk 



-irko 



ik — VFk\\ 



(3.24) 



Applying the residue theorem, the ko- integration yields a constant equal to —1 for k\\ > 
and zero for fcii < 0. Thus the A;- integration gives a result of order 1. The factor Aq - ^ is 
cancelled by the summation over sectors. Therefore the tadpole amplitude is O(g^). Direct 
computation shows that the turtle diagram is of the same order. Comparing with ( |3.16| ) we 
obtain 

Sp® = O(gW/\ ) (3.25) 

The renormalization of the chemical potential deforms the singular surface in the propagator 
- the Fermi surface — and modifies the original form ( ft .ICQ of the effective action. This 
will be cured by a change of variables discussed at the end of this section. 

Let us now turn to the renormalization of the couplings g. To determine g^\ we have to 
sum over all connected graphs with four external legs. In fact, the graphs that are not 1PI 
give a vanishing contribution. These graphs are of the form 




(3.26) 



As in the renormalization of the electron propagator, the amplitude of such a graph vanishes, 
as a consequence of momentum conservation. 

The tree level contribution to g^\uji, . . . ,#4) is just g(°'(ui, . . . ,£4). The one-loop cor- 
rection is described by the diagrams 




+ 




uj, k 



(3.27) 
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The first diagram corresponds to the sum of the following three graphs 




(3.28) 



The amplitude is of order at most 0((g^) 2 / Xq). In fact, for uji ^ ujf, the number of 
cJ's contributing to the ^-summation is 0(Aq -2 ), due to sector momentum conservation. 
Consequently the amplitude is of order 0((g^) 2 / Xq). But when uli = ulf the number of 
terms in the sum is actually O{X l ~ 1 ). We thus must take a closer look at the integration 
involved in the calculation of the amplitudes. The fc-integration extends over k's satisfying 
^ < \ulk\ < k F and < \u;(k + Pi — f>f)\ < k F . Explicitly the loop integral is given by 

/ d d k I dk - - 

J J ik -v F u)k i(k +p it0 -p ffi ) -v F uj{k + pi-p f ) 

In order for the integral to be different from zero, the two poles must not lie in the same 
(upper or lower) half plane, i.e., ulk and u(k + Pi — Pf) must have opposite signs. The set 
of fc's satisfying all three conditions is of measure zero since \uj(pi — pf)\ < Hence the 
integral vanishes. Thus the amplitude of the graph is of order 0(((7(o)) 2 /Aq), as claimed, 
except for exceptional configurations of external sector momenta (| uji — ujf |~ ^). 
Next, we turn to the second diagram in (|3.27|) corresponding to the graph 



< (3.29) 



UJ • . p 



When Ui ^ the number of nonvanishing terms in the cJ-summation is 0(A^ 2 ), and the 
amplitude is at most of order 0((g^) 2 / Xq). But, for the BCS configuration, uli = — cJ-, the 
situation is completely different. There are now 0(Aq _1 ) possible choices for the internal 
particle momenta, and no miracle makes the amplitude vanish. In fact, the loop integration 
for pi = pf = is 

d d k [ dk — (3.30) 

J kl + {v F ukf 

which is strictly positive. Thus the amplitude is of order 0((^^) 2 /Ao _1 ). 



B,7 
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Inserting the correct scale factors, we get the following flow equations for the quartic 
couplings: 

g W = g (P) + O((g^) 2 /X ) , for + -4 (3.31) 

9$x = & + 0({g^) 2 ) = (3.32) 

Thus, for the BCS couplings, the tree level does not yield the complete contribution to lowest 
order in -^-! In order to understand the flow of the BCS couplings, we have to investigate 
the loop corrections in more detail. This will be done in the next section. 

The iteration step from scale j to scale j + 1 is analyzed in a similar manner as from 
scale to scale I. Neglecting vertices of degree > A imp* and ip, we obtain 



z U+i) = + 0((g^) 2 /X j ) (3.33) 

= vjp + OUgVy/Xi) (3.34) 

= 0{g {j) /\ 3 ) (3.35) 

gii+D = g U) + Odg^/Xj) , for uj, ^ -uj[ (3.36) 

9^ = & + 0((^ s ) 2 ) , for % = -<% (3.37) 



We observe that the couplings g( j \uj 1: ...,o>4), for uoi ^ —a;-, do essentially not flow. But the 
BCS couplings may change considerably under the renormalization group flow. If the BCS 
channel is turned off, perturbation theory is valid, and the system approaches a Landau- 
Fermi liquid. If g^cs grows i n j then perturbation theory breaks down. In this situation 
the system becomes a superconductor, as studied in more detail in the next section and in 
Chapter 4. The unlimited growth of ^bcs i n 3 then reflects the fact that we are performing 
a perturbative analysis about a state that is not a ground state. As a matter of fact, 
superconductors do not possess a Fermi surface. 

At this point, a comment on contributions of degree > 4 in ip*, ip to the effective action 
is appropriate. The decimation of degrees of freedom (i.e., the integration over degrees 
of freedom) in an iteration step j, j = 0,1,2,..., of the renormalization group procedure 
produces terms like 

2k 



1 




(3.38) 



of scaling dimension — (d— l)(k — 1) in the effective action . The form and renormalization 
flow of the dominant contributions to these terms can be studied quite explicitly. In the next 
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iteration step, from j to j + 1, these terms yield contributions to the dimensionless coupling 
constants </ J+1 ) of the terms of degree 4 in S^ 1 ^ that turn out to be of order 1 = (t-)°, 
even for incoming sector momenta cJj 7^ — u^, 



2k 





(3.39) 



(jj; 



CO 



and hence may be important. They are, however, of higher order in g^\g^~ l \ .... A careful 
analysis (presented elsewhere) reveals that these contributions induce a finite flow of the 
couplings towards RPA-(random phase approximation) type fixed points. For short- 
range two-body interactions, and if g^ <C 1, this represents an unimportant modification of 
( 1-3.36 ). However, for long-range (e.g. Coulomb) two-body interactions, or if g^ is not small, 
the modification in the renormalization flow of the g^'s due to the terms of degree > 4 in 
if}*, ip in the effective action is essential and is intimately connected with the phenomenon of 
screening. A detailed discussion of these matters would go beyond the scope of these notes. 

Next, we analyze the renormalization of the chemical potential. At scale j, the effective 
action is given by (we use unsealed fields) 



S. 



(?) 

eff 



(T,UJ 

+ higher degree terms 



Su 



CO 



0) 



, , of the fey-variable. Note that ^-y 



k F 0( 



,0) • 



We may perform a shift, fen — > k 

Thus, if the coupling remains approximately constant under the RG flow and sufficiently 
small, the shift of the fen-variable will always be smaller than 4 s . The integration measure 
is invariant under this coordinate transformation. The transformed action reads 



qU) 



.^(i)-i^ J d d+l kipl a (k)(ikQ - v^k^ipz^k) + higher order terms 



Hence S$ is again of the form ( p.lOj ) (except that the domain over which fe is integrated 
has changed slightly). 

Our results motivate the following interpretation of the renormalization group trans- 
formation: Assuming that the interacting system has a Fermi surface (i.e., that the BCS 
channel is turned off), the RG procedure enables us to approach the renormalized Fermi 
surface stepwise, starting from the Fermi surface of the noninteracting system. In the limit 
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Figure 4: Renormalization of the Fermi surface 



j — > oo, we reach the physical Fermi surface of the interacting system. Imposing renormal- 
ization conditions on Z^°°\ Vp°\ amounts to solving a "final condition problem" in 
the space of running coupling constants. The RG recursion generates the discrete dynamics 
of coupling constants, as the scale parameter j is varied, but instead of initial conditions, 
experiment fixes final conditions. What we actually do is that we integrate the equations for 
the running couplings backwards from the Landau-Fermi liquid fixed point (j = oo) down to 
j = 0. The values /i^, Z^°\ Vp\ and parametrize a microscopic system, whereas 
Z(°°\ Vp and g^ parametrize a universality class of scaling limits describing macroscopic 
states. 

3.3 The BCS channel 

In the previous section we observed that the one-loop correction in the renormalization 
group flow of (?bcs is of order zero in i; gees might be a relevant coupling. To reach a better 
understanding of gscs, we determine the flow of ^bcs, taking into account all corrections of 
order zero in j (but ignoring terms of degree > 4 in ip* and ip in the effective actions; see 
( 13.381) , ( |3.39|) ). We omit the scale index when no confusion may arise. 

First, we have to determine all four- legged diagrams corresponding to corrections to <?bcs 
of order zero in ~. These are precisely those diagrams whose amplitude is of order O(j^j). 
Every graph with four or more external legs can be constructed from a uniquely determined 
diagram containing no two-legged subgraphs. This is done by replacing each inner line by 
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an appropriate graph with two external legs. We denote by 




(3.40) 



the sum of all graphs with two external legs; i.e., in a diagram we have to replace each such 
symbol by the sum of all two-legged diagrams. Let us consider an arbitrary four- or more 
legged diagram with no two-legged subdiagrams, with all inner lines carrying this modified 
electron propagator line. We shall compute the amplitude of such a graph by just replacing 
the electron propagator by the amplitude of the modified electron propagator line. (To see 
why this is correct, we expand each modified propagator into the corresponding sum of terms. 
Every summand is the amplitude of a specific diagram to which a combinatorial factor is 
assigned that is, in general, smaller than the combinatorial factor belonging to that diagram. 
But one easily proves that summing up all terms corresponding to a given diagram yields 
the correct factor.) 

What are the connected diagrams with four external legs whose amplitude is of order 
0(js=t)7 Consider a four- legged graph with n interaction squiggles. The squiggles provide a 
factor of order A"""^ -1 ). There are 2n— 2 inner lines. But the ^-functions of sector momentum 
conservation restrict the sum over the cJ's to n — 1 independent sector momenta. Thus, in 
order for the amplitude to be of order O(jj^j), all the n — 1 sector momentum summations 
must extend over the entire (d — l)-dimensional Fermi sphere! In this case, the diagrams 




(3.41) 



with n > 0, are of order O(jjrr), unless a miracle happens that makes some of them vanish. 
Using the above argument, one can prove that the graphs ( |3.41| ) are the only four-legged 
diagrams with no two-legged subdiagrams which are of order 0(t3=t)- 

We conclude that if we do not include terms of degree > 4 in if)*, ip in the effective action 
then all the diagrams that yield a correction to gBcs of order zero in j- are contained in the 
set of graphs built from ( |3.41| ) by replacing all the inner lines by ( |3.4U| ). In fact, in ( |3.40 ) 
only the graphs of order zero in ~ have to be taken into account. Let us consequently define 




(3.42) 



33 



as the sum of all two-legged graphs corresponding to corrections of order zero in j- to the 
quadratic part of the effective action. The diagrams that have to be computed are 



(3.43) 



with n e IN. In our computation we just have to replace the electron propagators in ( ft.41| ) by 
the renormalized ones in ( |3.42| ). It is well known that the renormalized electron propagator 
is given by 



amplitude of — »— Os, kj— *— = — (ifc — VFtok + amplitude of — (lPlV* — ) (3.44) 

where only amputated 1PI graphs with amplitude of order zero in j- are included on the 
right hand side. The results discussed in the last section then show that 

amplitude of — *—u3,k\— = — (ik — vpCuk + A5/ii) _1 (3.45) 



where 5/j.i is the renormalization of the chemical potential, and only contributions of first 
order in -r are taken into account; 5fi± depends on Z, Vp and g. 

Now we are able to compute the renormalized value of <?bcs- The amplitude of ( |3.43| ) 
with zero incoming and outgoing box momenta is given by 

1 \ n+1 

E (-l) n /5 n 5f B Cs(^,^n)fi'BCs(^n,^n-l) ■ ■ ■S'BCS^l,^') 



X d-1 

Wl,...£J n 

and j3 is a strictly positive number given by 

1 1 



(3 = J dkj_ dk\\ dk 



ik — Vpkn + \8jii —ik — Vpk\\ + A5/xi 



dkx dh kl + (v F k] - AW > ° (3 ' 46) 
We find that the renormalized value of c/bcs is 

OO / 1 \ n 

»Sjs(^^ + E t*=T E (-lr^s(^^)•••^s(^^ , ) (3-47) 



-' \ \ a 

=i V A j 

+0 



9 
A 



•j 
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In order to get a more explicit expression for the flow equation ( p. 47] ), it is useful to expand 



the BCS couplings, g BC s{^^') = 9bcs(^(^^'))i spherical harmonics 

oo 

g B cs(^, uj') = E 9ih{u, uj') (3.48) 
1=0 

with 

f nh , cqs(U(u},u')) d = 2 

hiU2,uri:={ £t '°\ \ (3-49) 

Pi, I G IN, are the Legendre polynomials. For d — 2, we assume 0(2)-invariance of the 
potentials (instead of only SO(2)-invariance). The normalizations of the functions hi have 
been chosen in such a way that 

da(tii) hi(u,uji)h v (uJi,Q') = 5i t i>hi(u,u') (3.50) 



|o)|=l 

Because 



TirE(-)= / /-(-)(■) (3.51) 



the r.h.s. of the flow equation (|3.47|) becomes (we omit the sub- and superscripts j) 

oo „ 

g BCS (uj,u) + J2{-l) n f3 n / d<r(ui) ■■■ da(uj n ) g B cs{^, <3 n ) '•• S'bcs^i, u) 

oo oo 

= 9M$, + E(- 1 )"/ ?n E 9l Q ■ ■ ■ 9iJio,h ■ ■ ■ Sin-i^Kfa 



n=l 

oo oo 

= E(E(- 1 )W +1 )^(^^ / ) 

/=0 n=0 

oo 

up to terms of order j-. The flow equation for the BCS couplings hence takes the form 

= r^-m + ° (f) ' /GIN (3 - 52) 

where /3j is positive and approximately independent of j; it depends on v^f 1 and 
The flow equations for different angular momenta, /, decouple to lowest order in j. But at 
order t, the flow equations for different Vs are coupled (Kohn-Luttinger effect). 

Before we analyze the flow of the couplings, we study the role played by electron spin. The 
potential v(x) was supposed to be independent of spin. Thus the couplings g^\uJi, ■ ■ ■ ,£4) 
are spin- in dependent, too. We thus divide the quartic term in the effective action into a 
spin-singlet and a spin-triplet part. We define 

s 1 
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Clearly <\>%^ a%a , = T^^a',^ and hence 4>h u ^,o,c = (0 s wil1 correspond to spin-singlet 
pairing, and 0* to spin-triplet pairing). The quartic term in the action becomes 

lx^ z ~ 2 E / d d+lk i ■ ■ ■ dd+lk ± (^) d+1 Kh + h-k 3 - k 4 ) 

<3 l+ l3 2= l3 3+ l3 4 
<J,cr' 

+^(g(ui,U2, u 3 , u*) - g(ui, u 2 , c3 4 , ^))^kA k ^kA k ^k,^,a'A k ^ fc *)} 

Thus, we define the singlet and the triplet couplings as 

s 1 

^(wi,^,^,^) := -(#(^i, £2,^3,^4) ±#(£1,^2,^4,^3)) (3.54) 

Because of the property that g(cJi, u3 2 , u5 3 , u; 4 ) = <?(o>2, ^4, ^3) we see that 

s s s 

5f*(cJi,cj2,^3,cj 4 ) = (£2,^1,^3,^4) = ±g\u 1 ,u 2 ,u 4: ,U3) (3.55) 
Using these symmetry properties, the quartic term in the effective action is found to be 

\^=i z ~ 2 E / dd+lk ^ ■ ■ ■ d<i+1 ^ (^) d+1 Kh + h-h- h) 

a, a' 

{g s (£i, u 2 , £3, ^Wala^a'A^ WskfrS A k *> k ±) ( 3 - 56 ) 
+{g t (u 1 ,u 2 , w 3 , W^^fei, h^k^a'Ah, h)} 

i.e., the quartic term is the sum of a singlet and a triplet part. The coupling of the singlet 
part is g s , and the coupling of the triplet part g % . In terms of the original BCS couplings, g s 
and g 1 are given by 

s 1 

#bcs(^> = -(^bcs(w, £') ± #bcs(^, -w')) (3.57) 

Using that Z(cD, — cD') = + 7r, the functions /ij can be seen to satisfy hi(uj,—uj') = 

(—l) l hi(£u,£u'). Consequently the singlet and triplet BCS couplings are 

00 

#bcs(^0 = Y,92qh 2q {u,u') (3.58) 

q=0 

00 

0BCs(^') = E^2g+lH+l(^^') ( 3 - 59 ) 

In the expansion of (7bcs> on ^y even an gular momenta appear, and the expansion of <?* BCS 
only involves odd angular momenta, as required by the Pauli principle. 
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We now return to the analysis of the renormalization group flow. The flow equations 
( |3.52D can be written in the form 

«r>-a w = -^4 + o(*H (3.60) 

Instead of studying this difference equation we propose to investigate the corresponding 
differential equation. The differential flow equation is obtained in the limit where the size, 
M — 1, of the scale change in an iteration step tends to zero. Let us write = g\ ■ , and 
similarly for all the other running coupling constants. We define 

9i{t) := gf Xo) , g(t):=g^ , Z(t) := , v F (t) := (3.61) 

Consider a scale A = e*A and set M = e*' - *, if > t. We divide both sides of ( |3.60| ) by if — t 
and take the limit t' \ t. The l.h.s. yields 4;gi{t)- The coefficient (3 = /3(t',t) vanishes in 
the limit if \ t, and only terms linear in t' — t have to be kept on the r.h.s. of fl3.60p . Thus 
we linearize fl(t',t) : 

(3{t\t) = {t'-t) 1 {t)^0{{t'-tf) (3.62) 

7 (t) = £pPW,t) \t'=t>0 (3.63) 

The positivity of 7 follows from the monotone growth of /3(t',t) in t'. Contributions to (3 
corresponding to diagrams with two or more loops are of order 0((t' — t) 2 ). Thus, for the 
calculation of 7, only one-loop diagrams have to be taken into account; in particular, we 
don't have to calculate corrections to the electron propagator as we did in the calculation of 
(3. Taking the limit if \ t, the r.h.s. of ([3.60 ), divided by t' — t becomes 

- igi {tf + 0{e-*g(t) 2 ) 

For the error term the same argument as above applies, and we conclude that only one-loop 
diagrams contribute. We thus obtain the flow equations 

j t 9i{t) = -19i{t? + 0{e- t g{tf) (3.64) 

where 7 = j(t, Z(t), vp(t),g(t)) is positive, independent of I, and approximately independent 
of t. Of course, we could have found these flow equations by just looking at all possible one- 
loop graphs, namely 



+ 
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The first one yields the first term on the right side in the flow equation and the second graph 
the error term. 

We propose to determine the flow described by fl3.64|) , neglecting the error term, and, in 



accordance with the results of the previous section, assuming that Vf, Z and g essentially 
do not flow. Then 7 may be approximated by a positive constant 70. Thus we have to solve 
the differential equation 

-j t 9i = ~logf (3.65) 
The solution of this equation is easily found to be 

»<*> = TTli ( 3 - 66 > 

Note that if gi{0) < the solution blows up at a scale 

1 1 

t = — , i.e. , for A = e^^A (3.67) 

7o0i(°) 

Thus if there is an angular momentum channel, I, with attractive interactions ((7/(0) < 0) 
the flow seems to diverge at a finite value of t. But this just means that perturbation theory 
breaks down when t ~ — (70^(0)) -1 , and we shall have to employ nonperturbative methods. 
The failure of perturbative methods is due to the circumstance that we are expanding around 
the wrong state ! 

Seemingly, everything is fine if gi{0) > V/. However, we have to remember that in 
( |3.65D we have omitted the error term which couples the flows of the running couplings g\ 



for different values of I. Without fine-tuning of the microscopic two-body potential, it will 
typically happen that gi(t) < 0, for some I, at some scale t. At that point, gi(t) will start 
to grow untill perturbation theory breaks down, and the ground state of the system will 
be superconducting. Thus, generically, a rotationally invariant system of non-relativistic 
interacting electrons will be a superconductor for small enough temperatures. This is the 
celebrated Kohn-Luttinger effect first studied in this fashion by Feldman, Magnen, Rivasseau 
and Trubowitz. 



4 Spontaneous breaking of gauge invariance, and su- 
perconductivity 

In the last chapter, we have argued that if some BCS coupling becomes negative at some 
scale A J0 then g\^ grows in j, for j > j , untill it becomes so large that our perturbative 
treatment breaks down. This phenomenon is the signal for an instability of the RG fixed 
point around which we are doing perturbation theory. Indeed, the Landau-Fermi liquid state 
is not the true ground state of the system anymore, and we expect that the RG flow drives 
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the system towards a new RG fixed point describing a superconductor. The global U(l)- 
symmetry of non-relativistic many-body theory (gauge invariance of the first kind) turns out 
to be spontaneously broken in the new (stable) ground state of the system. 

The purpose of this chapter is to analyze superconducting ground states and the asso- 
ciated breaking of gauge invariance. This is not an entirely simple story, and we therefore 
focus our attention on the simplest example, that of an s-wave (BCS) superconductor. 

Thus, we consider a system with the property that, at some scale Ao 3> 1 (with |<?(Aol 
1), gl% < and \g\ [ | <C —g^-, for / = 1,2,3,... . According to the results of Chapter 
3, there is then some j sc ps ^ , see ( |3.67| ), such that at scale Xj sc > A <7o <C 0, 

\g\ lsc) \ < -9o sc \ for I = 1,2,3,... , and (for (2 X ± -Cu 2 ) \g^(uj u . . . ,dJ 4 )| < -<7?' sc) - This 
suggests that we neglect all terms of degree > 4 (in ip*, if)) in the effective action of the 

Jjsc) 

system at scale Aj sc , except the s-wave BCS term with coupling constant ripr- The resulting 

jsc 

effective field theory is the one first considered by Nambu and Gorkov. 

There is some useful notation to be introduced: Let [uj] be the ray through the origin 
containing uj and — uj. We might think of uj and — uj as the two chiralities of a 1 + 1 dimensional 
system of relativistic fermions. We define field variables 

- ( & ) • - ( ) <4i) 

^[o3]T := (V'-^V'St) > ^[<3]i := (ip-au VtoJ ( 4 - 2 ) 

We thus group together field variables belonging to sectors on the same ray. Ordered ac- 
cording to spin indices, we consider them as entries of a four-component quasi-particle field 

^ : =(tZ) ' ^ := *t'W (4-3) 



named after Nambu and Gorkov. Let Vyi denote the two-dimensional, complex vector space 



whose elements are of the form ( ^ ) , and let V sp i n denote the two-dimensional, complex 



vector space whose elements are SU(2) spinors ( ^ j . We can think of the four-component 



object defined in ( f4.3| ) as being an element of V sp i n ® Vp]- Stressing analogies to 1+1 



dimensional, relativistic models ("dimensional reduction"), we also define gamma matrices 
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<7i , 7 1 :=l 2 ®a 2 (4.4) 



((7i = ( ^ q ] and o"2 = ( n * ] are Pauli matrices). Then we see that 



V>[fiJ]T = #]T°i » V>[<3]1 = ( 4 - 5 ) 
= ^7° (4.6) 
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The effective action at scale j sc , simplified by omitting all subleading terms of degree 2 
(curvature of Fermi surface) and > 4 (g^ sc \ g\>Q ~ <C —g\=o >■ • • ) is given by 

m 

+ TTd^T E / [^[fiJ^l ® l2)^p]^P']((Xl ® 12)^13'] 
4A i sc [<3],[<3']' / 

+ ^ (<r 2 ® l 2 )Vrf[<3'] (o- a ® l 2 )^p']] ( 4 - 7 ) 
It only includes the s-wave BCS interactions and has the following basic features: 

1) The action strongly resembles the one of a 1 + 1 dimensional, relativistic quantum 
field theory with N = N 3 * c) fermion flavours and quartic self-interaction, such as the 
chiral Gross-Neveu model. Its infrared properties are closely related to those of the 
chiral Gross-Neveu model. The perturbative infrared renormalization around the Fermi 
surface and the -^-expansion are virtually identical. 

2) The action, written as in (|4.7| ), exhibits a manifest global U(l)-symmetry, given by 
with a G [0,27r). 

The generator of this symmetry is the particle number operator N; it is the usual gauge 
symmetry of the first kind. In the superconducting phase of the system, this symmetry 
is spontaneously broken. We recall that the Mermin- Wagner theorem states that, for a 
field theory model in d + 1 = 2 dimensions with a continuous symmetry such as the chiral 
Gross-Neveu model, the continuous symmetry cannot be broken spontaneously. In the chiral 
Gross-Neveu model, continuous symmetry breaking only occurs in the N — > oo limit. In 
spite of the formal similarities between nonrelativistic many-body theory in d > 2 space 
dimensions, in the Nambu-Gorkov approximation, and the chiral Gross-Neveu model in one 
space dimension ("dimensional reduction"), the Mermin- Wagner theorem does actually not 
apply to the former, and spontaneous breaking of the U(l) gauge symmetry is possible in 
many-body theory (for d = 2 at temperature T = 0, and for d > 3 at sufficiently small 
temperatures) ! 

By direct calculation we see that 

i) [Sj] {o- x ® i 2 )V>[o5] = {V'-arV'fiu + ^h^-zi + ^-oifai + ^ai^-an) 

= 2{0i_^ iT (x, x) + ^, u (x, x)} (4.8) 
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and 



= 2z{^_ £3>iiT (x, x) - *)} (4.9) 

where s has been defined in ( p.53| ). Thus, these quadratic expressions correspond to real 
and imaginary part of the BCS order parameter, respectively. 

The action Sesty,ip) defined in ( |4.7| ) can be replaced by a more convenient, equivalent 
action, S, that is quadratic in -0 and 0? and depends on a complex Lagrange multiplier field 
= 0i + ^02- The new action is given by 



Sty, ^,0,0) = E / (7°^ - v Fl l uV)^ [a] \d d+1 x 
[<3] J 

+#E / fe^l ® 1 2)^M01 - ^[u5](^2 ® 1 2 )^M0 2 ]^ +1 X (4.10) 

+\ J (<1>1 + <Pl)d d+1 x 



Use) 



where 2g 2 = — ^ d _ 1 > 0. We note that, under the U(l) symmetry discussed in remark 2) 



JSC 



above, the field transforms as — > e 2ia (p. We emphasize that S e sty,ip) and Sty, ip, 4>, 0) 
are equivalent in terms of their physical content: It is easily checked that, after functionally 
integrating out the 0-field, 



J VfaVfae-^'*'*'® = const e 



the original action S e ^ty, ip) is restored in the exponent. But Sty, ip, 0, 0) is much easier to 
work with, because it is quadratic in the fields ip and ip. The Bose field will turn out to 
describe the Cooper pairs of electrons. 

The interaction between fermion- and Lagrange multiplier fields is described by 

9^2 [0(a:)^],T^M,| + ^{x)^ [s]> ty [a]A } d d+1 x (4.11) 



This expression shows how interaction vertices between electrons (holes) and bosons 0, 
are organized in the Nambu-Gorkov theory. One Nambu-Gorkov vertex is equivalent to the 
following four vertices (each arrow points to a ipa,<r)'- 



,0 
i 



g 



-g, j^g , | 



(4.12) 
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One can join only the first two vertices with each other and the second two vertices with 
each other. There are two possibilities of joining the first two vertices : 



10 10 10 10 

^'J + "^"L 

g g g " g 

and two possibilities of joining the second two vertices, obtained by exchanging j and | . 
Along a string of electron (hole) propagator lines, the different sector labels Cu, —u and boson 
fields 0, have to occur in an alternating pattern. Because the first two and the second 
two vertices never intercombine, the spins can be omitted if a factor of 2 is attached to each 
loop. 

In the following, we focus our attention on the boson fields, and we will attempt to 
eliminate the Nambu-Gorkov fermions. Because S(ip, ifj, 0, 0) is quadratic in -0 and ip, one 
can perform the fermionic functional integration explicitly: 



J V^V^ e~~ s ^^ M) = exp (-^ J d d+1 x |0(x) 



■ det[( 7 °«9, - vrfuti) + g ( £ x) *W ) ® 1 2 ] (4.14) 

After normalizing the right hand side by division through det[y dt — fi?7 1 u;V], the determi- 
nant 

can be evaluated using the identity 

det(l + A) = exp[]T ^ ± Tr(A n )] (4.16) 

n=l n 

The term Tr (A n ) is the amplitude of the n-th order one-loop diagram 



(4.17) 



42 



/ o (b(x) \ 

with a factor of {g I ^, . J ® ^ 011 eac ' 1 ex ^ ema l li ne an d a Nambu-Gorkov propa- 
gator on each segment of the loop. From our discussion of the possible pairings of Nambu- 
Gorkov vertices we conclude that loops with an odd number of vertices vanish. The expres- 
sion in the exponent on the right side of ( |4.16|) reduces to — ~ J2%Li \Tr(A 2n ). Calculating 
Tr (A 2n ) in terms of the original electron (hole) propagators, the spin summation can be 
absorbed in an overall factor of 2 cancelling the | in the exponent. For a fixed [u], a 2n-loop 
looks like 

(4.18) 




inserting alternatingly electron (hole) sector labels uj, —uj and boson fields <fi and along the 
loop line. In order to compute the full amplitude of such a graph in momentum space, we 
have to integrate over the loop momentum (ko, k) for arbitrary external momenta (pj,o,Pj), 
j = 1, 2n, which is a difficult task. 

It will turn out to be useful to expand the amplitude corresponding to ( |4.18j ) into a 



sum of terms that look somewhat more manageable: Anticipating spontaneous symmetry 
breaking, we assume that, in a (superconducting, extremal) ground state of the system, the 
Bose field has a non-zero expectation value C . The modulus of |0 C | is determined by the 
values of physical parameters (the density of the system, the strength of g$ , etc.), while 
the phase of C , an angle in [0, 2n), is only fixed after suitable symmetry breaking boundary 
conditions have been imposed — as usual in the study of systems with spontaneously broken 
symmetries. We thus decompose the Bose field 4> into a constant part, <fi c , and a fluctuation 
field, x(x): 

(f)(x) = <fi c + x(x) 

${x) = C + X(x) (4.19) 

The field xi x ) describes small fluctuations of the Cooper-pair condensate around <j) c . The 
decomposition of 0(x) induces a decomposition of the amplitude corresponding to (|4 . 1 8|) into 



a sum of monomials in x an d X- F° r each fixed n and [u], this decomposition, a binomial 
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series in x(p) an d xip)i can be described pictorially, as follows: 




+ terms of degree > 3 in x an d X- 



Each dot without an external leg stands for a factor g<j) c (2ir) d+1 5(p) or for its complex 
conjugate. The external momenta flowing into the diagram at all such dots vanish, because 
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(j) c is constant in x-space. The graphs with one external leg are summed over all 2n possible 
positions of either %(0) or x(0) on the loop. The corresponding amplitudes (which are linear 
in x or x) vanish if <ft c is the expectation value of <fi in the ground state of the system. The 
effective action for the field x(j?) is thus given by a power series in x an d X '■ 

S(X, X) = 5 (0) + X) + S {2) (x, X) + - (4.20) 

where S^ r \ r = 0, 1, 2, is a sum of monomials of degree r in x and x- We remark that the 
first three terms on the right side of ( f4.2U ) also contain the contribution 



+<P c X(p)(2Ti) d+1 5( P ) + ^ cX (p)(2n) d+1 5(p) + x{p)x(p)) (4-21) 



from the part of the original action S ( defined in ( |4 . 1 0| ) ) quadratic in and 0. 

After dividing by the total volume of the system, i.e., by (2n) d+1 5(0), the amplitude of 
the 2n-th order loop without any external x~ an d X-^ e S s , an d f° r a fixed ray [u], 




(4.22) 



is given by 

2 f d d+1 k( ~ 9 ^\ y (4.23) 
JrxBj kl + v 2 F (Qk) 2 

Summing over all orders 2n ( see ( t4.16|) ) and all rays [u] , our result for is found to be 
S® = i f . d d+1 k\^ c \ 2 ((2n)6(k)) 2 

+2 W _ <t*»k £ - l2 ~ 9 Z e L 2 (2«)*»6(p = 0) (4.24) 



(2n)^5(p = 0) ( 1|0 C | 2 - 2^ / . d d+1 k log (l ' " 



2|J> 12 



The effective potential, U c g((j) c ,(j) c ), is defined as the density of S^, i.e., as S^ ) divided by 
the total volume of the system. Thus 

U^ c ) = \W - 2E L n log (l + , 2 f'f C L 2 ) (4-25) 

2 M ^Rxfla V k 2 + V 2 F {uk) 2 J 
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In the analysis of spontaneous symmetry breaking, the effective potential plays an im- 
portant role. In the approximation of mean field theory, the expectation value of the field <fi 
in an arbitrary, extremal ground state of the system is given by a minimum of U e s(4> c , <f) c ). 
Thanks to the minus sign in front of the integral on the right side of ( |4.25| ), the graph of 
the effective potential has the shape of a Mexican hat. The minima of U e R are obtained by 



setting the derivative of U e g with respect to 
values of g 2 ^ 1 ) 

0c ~ exp( 

g 



to zero. The result is that (for small 



TTVp 



g 2 (k F X J 



(4.26) 



where v F is dimensionless (in our units), k F has a dimension of inverse length, and g has a 
dimension of (length)", so that the dimension of \(f> c \ is that of (length) ~ , as it should 
be in view of the last term on the right side of ( |4.10| ). 



We recall that 2g 2 



(jsc) 



(see ( |4.10| )) and hence 
\^2k F v F \, 



d-l 
2 

Jsc 



— n la 

go 



cxp( 



2tiv f 



g { o 3s 



>Ud-l 



(4.27) 



\9$' 



=)i 



(kp_ 



d-l 



for small values of 

At values of C minimizing U e s( 
= 0. (The equations -MH 



i.e. 



is given by ( |4.26| ).) From now on, 



, C ), the terms in S(x, x) linear in x or x must vanish, 
- and = are, of course, equivalent; the solution 
will denote a minimum of U e s- 

and 

XX> respectively. They can be found by calculating the amplitudes corresponding to the 
following sums of diagrams: 



The term S^x, x) consists of three different contributions, proportional to \ L \ X 2 



E 



all positions 
of two external x _ie g s 




x(-p) 
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We are interested in calculating the Taylor series expansions of the amplitudes corresponding 
to these sums of diagrams in p around p = 0. The terms constant in p yield the coefficients 
of the "mass terms" of the fluctuation field \. These coefficients are found to be given by 

tfFcY i d d+lk 1 (4-28) 

(coefficient of term proportional to x 2 )i the complex conjugate of this expression (coefficient 
of term proportional to x 2 ), and by 

1 - 2, 2 V [ d d+1 k k l±30L 

(coefficient of term proportional to XX) i where the "|" comes from the term | J d d+1 x XX m 
the original action S (see ( |4.10| )). 

It is convenient to introduce polar coordinates in field space by setting 

cf>(x) = \<p c \e w + (Xt(x) + iXi{x)y 6 (4-30) 

where the phase 6 of the ground state expectation value of 0(x), ~ <fi c , is fixed by the 
boundary conditions imposed on the system. The component Xt of x is parallel to <f) c , i.e., 
it is transversal to the manifold of minima of U e s, a circle, at C . The component xi is 
perpendicular to C , hence tangential to the manifold of minima of U e s at 4> c . 

A straightforward calculation (involving ( |4.26| ), i.e., the fact that <ft c minimizes U e s) now 
shows that the mass terms for the field x combine to 

3/2 'd^xxl (4-31) 



2 

where 

M 2 = 16s 4 |0 c | 2 y / d d+1 k 1 (4.32) 

In accordance with the Goldstone theorem, the effective field theory for x contains a massless 
field xi describing the Goldstone bosons and a seemingly massive field xt of mass M, where 
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M is equal to the square root of the curvature of £7 e ff in the radial direction (parallel to C ) 
at <j) c . Actually, it will turn out that the degree 3 terms, of the effective action of x 
couple Xt to the second power of xi- Thus the field quanta of Xt are resonances that decay 
into pairs of Goldstone bosons. 

We note that all terms in the effective action for the field x that are local and do not 
contain any derivatives of x an d X can be determined from U e s((f) c + x, C + x), with x( x ) = 
const, by expanding in powers of x an d X- 

We can also compute the kinetic terms of the effective action of x- Let us denote their 
inverse propagators by U XtX (p), U XtX (p), U XtX (p) and IT XiX (p), where p is the external mo- 
mentum, and the indices stand for the external legs. For example, the coefficients, ct XiX and 
(3 XiX , °f the contribution / d? +1 p x(p)( a x,xPo + Px,x\p\ 2 )x(p) t° the effective action are given 
by' 

1 d 2 

a x,x = 2 ~dpj \p=° U x,x(p) ( 4 - 33 ) 
1 d 2 

Px,x = ^ Ip=o n x , x (p) ( 434 ) 

Pi {i = 1, .., d) is the i-th component of p. Due to rotational symmetry, (3 XiX is independent 
of the choice of i. All other coefficients can be expressed in terms of derivatives of IT XiX and 
n xx . In d = 2 space dimensions and for small values of g 2 kF, we obtain the following results: 

Q « - a ^ii^ (4 - 35) 

k A 

a xx = a xx ^ - F 3 ff ,- (4.37) 
x,x x ' x 24ttw f |0 c | 2 v 7 

^ = £ (4.38) 
x,x x,x 96tt|0 c | 2 v 1 

Splitting x into a transversal and a tangential mode, we arrive at the expression 

for the kinetic part of the effective action. For g 2 kp small and d = 2, the contributions to 
(x, x) °f order zero in the momenta turn out to be 



/ 



d 3 pid 3 p 2 4 kF T7 V ( - Xi (Pi ) ^ ( ~P2 ) X/ (P2 - Pi) 

3! 7T|0 c | d 

+Xt(Pi)Xt(-p 2 )x t (P2-Pi)) (4.40) 
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Rescaling Xt by a factor ( j^^fu p ) 5 an d X/ by a factor (4^7^2)2 , the lowest order 
terms in the effective action for the x-field take the following standard form : 

S(x,X) = J ^PXt{p){pl + jV F \p\ 2 )xt(-p) 
+ / ^PXi{p){pl + \vl\p\ 2 )xi{~P) 

+ J d 3 p 12v 2 F k 2 F X t(p)Xt(-p) (4.41) 

-»(P (1, )X«(-P (2) )*(P <2) ~P m ) + 3x t (p m )x,(-P <ly )X,(p m - P m 

At tree level, the propagator (xt(fyxt(x)) of the Xrheld in x-space decays exponentially 
in \x\, with decay rate M. Using expression ( ^4.41| ) to calculate radiative corrections, the 
behaviour of the propagator of Xt changes drastically. This is due to the vertex in S(x, x) 
proportional to x 2 Xt- The dominant one-loop radiative correction to the propagator of Xt is 
proportional to 

"d d+1 yd d+1 z ( X t(0)xt(y))o (xi{y)Xi{z))l (Xt(z)xt(x))o (4.42) 

where ((-))o indicates that the expectation value is calculated at tree level. Because the 
propagator (xt(0)xt(%))o decays exponentially in \x\, while (xi(^)Xi( x ))o ~ |a;| 1_d , expression 
( [4.42| ) is proportional to 

(Xi(0)xi(x))o ~| x \ 2 ~ 2d (4.43) 

Thus the seemingly massive field quanta of the field Xt are resonances that decay into pairs 
of massless Goldstone bosons, as announced. 

The form Q4.41 ) of the effective action of the x-field shows that in one space dimen- 



sion the fluctuations of xi are logarithmically divergent (logarithmic infrared divergence of 
/ <? 2 P p 2 +CO n St p i )- The prediction of mean field theory that the continuous U(l) gauge sym- 
metry is spontaneously broken is therefore wrong (Mermin- Wagner theorem). The same 
conclusion is reached in two space dimensions at positive temperature, T. The propagator 
of Xi m momentum space is then proportional to 

(4.44) 



(ksTn) 2 + const (p 2 + p?,) 



where the Matsubara frequencies n are integers. The term corresponding to n = yields 
logarithmically divergent fluctuations of Xh an d, again, the U(l) symmetry is restored. 
However, these systems exhibit a Kosterlitz-Thouless transition, as g is varied. 

But, for d = 2 and T = 0, or for d > 3 and at sufficiently small temperatures, the 
predictions of mean field theory concerning spontaneous symmetry breaking are qualitatively 
correct ! 
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The form ( [4.10 ) of the effective action S(ip,if},(j),(j)) and Eqs. (|4.8|) and (O) show that 
the field quanta of the x-field have electric charge ±2e, (where e is the elementary electric 
charge). If the Coulomb interactions between charged quasi-particles are incorporated in our 
theory then the Goldstone bosons acquire a positive mass ( "Anders on- Higgs mechanism"), 
as briefly discussed in Chapter 5 of Part I. 

The prediction of our theory for the mass, ~ g\<p c \ with \<p c \ given by ( |4.26| ), of quasi- 
particles of charge ±e agrees with the solution of the standard BCS gap equation. 
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